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We study the influence of the preparation of an open quantum system on its reduced time evolu- 
tion. In contrast to the frequently considered case of an initial preparation where the total density 
matrix factorizes into a product of a system density matrix and a bath density matrix the time 
evolution generally is no longer governed by a linear map nor is this map afflne. Put differently, the 
evolution is truly nonlinear and cannot be cast into the form of a linear map plus a term that is 
independent of the initial density matrix of the open quantum system. As a consequence, the inho- 
mogeneity that emerges in formally exact generalized master equations is in fact a nonlinear term 
that vanishes for a factorizing initial state. The general results are elucidated with the example of 
two interacting spins prepared at thermal equilibrium with one spin subjected to an external field. 
The second spin represents the environment. The fleld allows the preparation of mixed density ma- 
trices of the first spin that can be represented as a convex combination of two limiting pure states, 
i.e. the preparable reduced density matrices make up a convex set. Moreover, the map from these 
reduced density matrices onto the corresponding density matrices of the total system is afflne only 
for vanishing coupling between the spins. In general, the set of the accessible total density matrices 
is nonconvex. 

PACS numbers: 03.65.Yz, 05.30.Ch, 02.50.-r 



I. INTRODUCTION 



Closed systems are known to be an idealization. In general, real systems interact with their environment and exhibit 
properties that cannot be observed in finite closed systems, such as irreversibility of the time evolution and, related, 
the relaxation of observables toward stationary values and dephasing, or decoherence. Various techniques have been 
developed to treat the dynamics of open systems without explicitly considering the full Hamiltonian dynamics . For 
example, effective equations for the reduced density matrix of the considered open system, known as master equations, 
have been proposed long ago |^ H, |^ and still are of considerable interest because of their conceptional simplicity 
and potential usefulness [MM; New challenges in this field of fundamental physics have come from nanotechnology 
and quantum computing [g. 

Any equation determining the time evolution of a density matrix has to obey several general properties which 
guarantee that the density matrix stays selfadjoint, positive and normalized in the course of time. These general 
properties still leave much freedom and, in order to further restrict possible dynamical laws, additional requirements 
for the dynamics have been postulated 'O"! . One seemingly natural property that often is assumed without even being 
mentioned is the linearity of the time evolution, which generally is understood as a consequence of the linearity of 
the Schrodinger and the Liouville-von Neumann equation for closed systems. This argument, however, only works 
by analogy and no proof of the necessity of this requirement is available. Just on the contrary Pechukas [l^, [ill has 
recently shown that linearity may only hold if the initial state of the total system factorizes into a product of a density 
matrix for the open system and another one for the environment, and if a sufficient number of pure states can be 
prepared. The assumption of linearity is a prerequisite of a Markovian dynamics and of complete positivity [T^ Il3| . 
These properties then lead to the mathematically well characterized class of Lindblad master equations. From the 
physical point of view, however, these equations suffer from certain deficiencies 10]. They are restricted to the regime 
of weak coupling between the considered system and its environment. In particular the weak coupling assumption 
will fail at sufficiently low temperatures [13. [l5| . Moreover, there are various general statistical mechanical properties 
that are incompatible with the assumption of a Markovian dynamics [T^ . 

A few microscopic models of systems interacting with their environment can be reduced exactly to Lindblad mas- 
ter equations with time dependent coefficients [TtI I18| , thereby describing the single time non-Markovian reduced 
dynamics. 

The dynamics of an open system is determined by both, the full dynamics of the considered system interacting 
with its environment and the initial state of the complete system. The significance of the initial state was emphasized 
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in several works In an experiment this initial state is imposed by a preparation 

procedure. Here we will only be concerned with equilibrium preparations that lead to a thermal equilibrium of the 
total system in the presence of external fields that only act on the system and that are switched off finally. In this 
way, the initial state is described by a canonical density matrix: 



/■ = exp <[ -/3 I if - V F,X, I V , (1) 



where /? is the inverse temperature, H the Hamiltonian governing the dynamics of the total system, F ~ (Fj) 
are external, i.e. classical, fields, Xj the corresponding conjugate operators of the open system and Z = 

Trexp |— /3 ^iJ + FjXj^ I is the partition function of the total system. Here, Tr denotes the trace over the 

total system. In this way, initial states of the total system are reproducibly prepared. They correspond to a thermal 
equilibrium of the environment at a given temperature and to a state of the system that depends on the strengths 
of the external fields Fj. The set of density matrices that can be obtained upon variation of the fields forms the 
equilibrium preparation class. The reduced states belonging to this preparation class are determined by the trace 
over the environment, which is denoted by Tr^: 

pf-TrsP^. (2) 

The calculation of this trace is nontrivial in most cases and in general does not lead to the canonical distribution 
of the uncoupled system, ^ exp {— /3i/s} J_4, 28, 29, 30J where Hs is the Hamiltonian of the system in presence 
of the external fields and Zs the respective partition function. This particular form is only obtained in the limit 
of weak coupling between the system and the environment 0, l29l l3ll |. In the weak coupling limit, the equilibrium 
density matrix of the total system factorizes into a product of a system and an environment density matrix. This is 
an example of the factorizing preparation which leads to a product of a particular density matrix of the environment 
Pb and an arbitrary density matrix ps of the system: 

P'"' ^PSPB. (3) 

The factorizing preparation is assumed in most theoretical investigations though it is often difficult, if not impossible, 
to realize it experimentally. 

A mo re g eneral class of preparations has been suggested in the context of the path integral approach to open 
systems [2^: 

p^ ^Y^O^p^O'^, (4) 
i 

where p^ is defined as in the equilibrium preparation, eq. and Oj and O'j are operators that only act on the 
system's Hilbert space. For applications of this preparation class we refer the reader to Ref. |2^ . 

The state of the open system after the preparation results as the trace over the environment of the density matrix 
of the full system at that time, i.e. 

psW = TrBC/(t)pC/t(i), (5) 

where 

C/(i) = cxp I --^iJtj (6) 

is the unitary time evolution operator of the full system and p{0) = p the density matrix resulting from the preparation. 

Requiring an affine time evolution of the reduced density matrix ps {t) means that ps (t) is the sum of a term linear 
in, and one independent of ps, i.e. 

Ps(t)^T{t)ps + I{t), (7) 

where T{t) is a linear operator and I(t) is independent of ps- According to eq. lO the reduced density matrix at time 
i is a linear image of the initial full density matrix p under the successive action of the unitary time evolution of the 
full system and the operation of the trace. In order to obtain a map from the initial reduced density matrix ps to its 
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value at a later time t we introduce the blow-up map R{ps) that assigns to each reduced initial density matrix ps one 
belonging to the total system: 

P = Rips). (8) 

Its particular form depends on the initial preparation of the system. Expressing the initial total density matrix p with 
the help of the blow-up map by the reduced initial state one obtains the reduced time evolution of the system: 

ps{t)^TrBUit)Rips)U^{t). (9) 

An afhne time evolution as described in eq. Q will result only from an afline blow-up map. So one may ask under 
which conditions the blow-up map is affine. 

Because a density matrix is a positive normalized operator, we require that the blow-up map R acts on a convex set of 
reduced density matrices, i.e. a set that contains with each pair psi and ps2 all convex combinations Xpsi + (1 — X)ps2 
for all < A < 1. We assume next that the considered preparation provides such a convex set of reduced density 
matrices. For any particular preparation one has to check this property. We then further may ask under which 
conditions R preserves the convexity condition, put differently under which conditions a full density matrix that 
corresponds to a convex combination of reduced density matrices again is given by a convex combination. Then the 
preparation class also forms a convex set. For such a convex blow-up map one obtains: 

R{Xps,i + (1 - A)ps,2) = Ai?(p5,i) + (1 - A)i?(/3s,2). (10) 
This then implies that the blow-up map is afhne: 

R{ps)^Lps + x, (11) 

where L is a linear operator that maps reduced density matrices onto full density matrices and x is an operator of 
the full system. A proof of this statement is given in the Appendix 

We note that the mixing parameter A has to be identical on the left and on the right hand side of the eq. (|10|l . This 
is a consequence of the fact that the trace of R{ps) over the Hilbert space of the environment must coincide with the 
system density matrix ps'. 

Trs i?(ps) -PS- (12) 

For the factorizing preparation (jS)) the blow-up map R is always linear. It simply acts as the multiplication by the 
reference environment density operator ps- In the case of a classical system dynamics the role of density matrices is 
taken over by probability densities defined on the respective phase space. Then, any preparation can be characterized 
by a conditional probability density p{x\xs) for the state x of the total system given the state xs of the system 
|l9l |2(1 l2ll | . The corresponding blow-up map R is then given by the multiplication with this conditional probabilit 
and, hence, is always linear. No such simple construction scheme is available in quantum mechanics and Pechukas 1 
showed that the factorizing preparation is the only one for which the blow-up map is linear, provided that a sufficient 
number of pure states of the reduced system are contained in the preparation class. For the convenience of the reader 
the precise formulation of the theorem and a proof is given in the Appendix IbI 

In the present work, we will consider the influence of the preparation of an open quantum system on its reduced 
time evolution by means of the simple example of two interacting spins. One of those is considered as the system, the 
other one plays the role of the environment. The second spin is only a very crude caricature of a true environment 
which clearly fails to cause dissipation or dephasing in the system because of its finiteness. Nevertheless, it suffices to 
illustrate the influence of the preparation on the subsequent dynamics of the system. 

We assume that the total system starts from an equilibrium preparation, i.e. that its initial state is described by a 
density matrix of the form of eq. (^). It will be shown that in general this preparation renders the time evolution of 
the reduced system nonlinear. 

II. TWO SPINS 

Both interacting spins <Ti = {erf, al,af) and (T2 = (af , erf, erf) with Pauli spin matrices a^^, a = 1,2, i = x,y, z, are 
of total length s = 1/2. The first spin, ai is considered as the system and the second one, a2 takes over the role of 
the environment. Every density matrix of the total system then assumes the form 

/3 = ^ (1 -I- Si ■ cTi -I- S2 • <T2 -I- cTi • C ■ 0-2) , (13) 
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where 

S„ = (cr„) (14) 

denotes the Bloch vector of the spin a — 1,2 and the matrix 

C = (cTi^Ta) (15) 

denotes the correlation matrix of the two spins. The dot-product denotes the scalar product in three dimensions, e.g. 
Si • (Ti = Sfaf + S\a\ + Sfaf. The reduced density matrix of the system is given by the trace over the environment 
(i.e. the second spin), and hence becomes 

PS = Tr2 p = ^ (1 + Si • (Ti) . (16) 

Here we want to study the opposite direction, that is to go from ps to p. In particular, we look for conditions 
under which the respective blow-up map R{ps) = p is convex. We recall that the blow-up map is determined by the 
preparation process of the system. In the process of a preparation the state of the system is controlled by external 
fields F that ideally act only on the system, as given in eq. ^ for the equilibrium preparation. For the considered 
spin one can think of static magnetic fields. For interacting spins the two Bloch vectors and the correlation matrix 
will depend on the applied magnetic field. 

We next investigate the requirement of the convexity of the blow-up map R. A necessary condition for this property 
to hold is the convexity of the domain of definition of R. This implies that for any pair of reduced density matrices 
p^*, fc = 1,2 that result from two different values of the field, all convex linear combinations can be prepared by 
means of another value F3 of the field: 



^XpI^ + {1-X)pI\ (17) 



with F3 being a uniquely defined function of the fields Fi, F2 and A. Using the general form of the reduced density 
matrix in Hlt)|) one finds a respective equation for the Bloch vectors, reading 

Si(F3) - ASi(Fi) + (l- A)Si(F2). (18) 

This means that the Bloch vector of the system must be a uniquely invertible function of the external field. For the 
equilibrium preparation this is the case because then the derivatives of the Bloch vector components with respect 
to the field components coincide with the correlations of fluctuations of the first spin. These derivatives form the 
elements of the susceptibility matrix which is an invertible matrix in thermal equilibrium. Hence, Si(F) has a uniquely 
defined inverse F(Si). We investigate the consequences of the convexity condition pO|l of the blow-up map for the 
spin system. Using eqs. (I17II and H13I) for the total density matrix of the full system one finds analogous relations 
from eq. (|10|) for both the Bloch vector of the second spin and for the correlation matrix with the same field F3 that 
results from eq. (fT5)l . i.e., 

S2(F3) = AS2(Fi) + (l- A)S2(F2), 

C(F3) = AC(Fi) + (l-A)C(F2). (19) 

These are non-trivial conditions which in general will not be satisfied. Expressing next the fields F^ by the set of 
Bloch vectors Si = Si(F,;) that result for the respective fields we find by use of eq. H18|l the following relations for 
the Bloch vectors of the second spin and the correlation matrix: 

S2[AS} + (1-A)S2] = AS2[S}] + (1-A)S2[S2], 

C[AS} + (1-A)S2] = AC[S}] + (1 - A)C[S?]. (20) 

Hereby we introduced the notation S2[Si] = S2(F) and C[Si] = C(F). From these equations it follows that both 
S2[Si] and C[Si] are affine functions, see Appendix^ Therefore they can be represented as 

S2[Si] = A-Si+B, 

C[Si] = D-Si+E, (21) 



where B is a constant vector, A and E are constant, second order tensors and D is a constant third order tensor. 
Hence, these quantities must neither depend on the applied field F nor on the Bloch vector Si. 
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III. AN EXPLICITE ILLUSTRATION 



We consider the equilibrium preparation for the following simple two-spin Hamiltonian as an example, 

H = -F,al + eal + 5(t?4, (22) 

Here we only allow for a field in the z-direction. The two spins interact by their i-components. We study the 
equilibrium preparation class at the fixed inverse temperature (3 that results if the field assumes all possible values, 

— cx) < Fz < oo: 

= Z-iexp{-/3i?}. (23) 

Because the Hamiltonian H commutes with the operator ct| one can diagonalize H in the eigenspaces of erf erf. This 
then yields the four eigenvalues £i 



£3 = -^{F.+ef+g^ £i = ^{F^ + ef+g^ (24) 
and the corresponding eigenprojection operators P^: 

1 / F — e a 



P, = - (^1 + alal - -^{^l + ^l) + jyi^l - ) , * = 1, 2, (25) 

^ ( ^ -Z _z Fz ^ / _z _z\ t 9 t _x _x I y y " 



= a{al ~ - ^2) + + <^M)) , * = 3, 4, (26) 

such that H = £iPi holds. The canonical density matrix p^' at the inverse temperature /3 is a mixture of the pure 
states Pi with the Boltzmann weights pi = exp{— /3fi} / [2 (cosh(/3fi) + cosh{(3£2))], i.e., 

P^^ = Y^P^P^ " 4 + ^i^^i + ^2.^72^ + C,.'!?^! + Cyyafa^ + C^afal) . (27) 

i=l 

The a:- and y-components of the two Bloch vectors vanish. The non-vanishing z-components read 

Si. = l3{F,T+{p£i,l3£s)-eJ^_{p£i,p£3)), 

S2z = P{F,J'-iP£i,P£3)-eJ^+{p£i,(3£3)), (28) 
where the auxiliary functions J-±{x,y) are defined by: 

^ , . ysinh(x) ± xsinh(u) . „. 

^±{x,y)= \ \ — TT^. (29) 

xy (cosh(a;) + cosh(y)) 
Finally, the non- vanishing elements of the correlation matrix C are given by: 

Cxx = ~^gT+{^£i,^£^), 
Cyy = (3gT^{(3£i,f5£^), 

^ cosh(/3gi) - cosh(/3g3) 
cosh(/3fi)+cosh(/3£3)' 

As already stated above, in the present case of the equilibrium preparation the z-component of the Bloch vector 
of the first spin is a uniquely invertible function of the field magnetic F.. This can be shown by inspection from 
eq. (1^ . see also Fig. We note that the equilibrium preparation contains the pure system states ps — \ {l ^ cr.) 
asymptotically in the infinite field limit F. — > ±00. If one also took into account fields that couple to the other spin 
components ax, Cy, the eigenstates of these components could also be prepared asymptotically. In the case studied 
here, however, no other pure states than the eigenstates of cr^ can be prepared. Thus, for this particular case one of 
the conditions of the Pechukas theorem to hold are not met. 
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FIG. 1: The Bloch vector component Szi as a function of the field PF^ resulting from eq. 12811 for /3e = 1. With increasing 
coupling parameter l3g = 0.5, 1, 1.5 the slope at K: = decreases. Apparently, the functions Szi{F) are monotonic and hence 
possess a unique inverse. 
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FIG. 2: The Bloch vector component Sz2 as a function of Szi resulting from eq. I|28|l for /9e = 1 and different values I3g = 
0, 0.5, 1, 1.5. Note that a strictly linear dependence results only for the case of vanishing coupling g = 0. An approximate 
linear regime exists in a neighborhood of Siz = 0. 



A. Testing convexity of the blow-up map 

We now come to the discussion of the eqs. H21f) which are necessary conditions for the blow-up map to be convex. 
Because in the considered preparation class only the z-component of the field is varied and because of the symmetries 
of the considered Hamiltonian these equations need only be checked as functions of 5*21 for the Bloch vector component 
and the diagonal elements of the correlation matrix, respectively: 



S^Szi] = ASzi + B, 

Cxx[Szl] — DxxSzl Exx, 

Cyy[Szl\ — DyySzl+Eyy, 

Czz[Szi] — DzzSzi+Ezz- (31) 

If the blow-up map were convex these equations would have to result from the eqs. (|28II3U|) by eliminating the external 
field Fz . There is no need to perform the lengthy calculation to see that the above equations hold only if the system- 
environment interaction is absent, i.e. if the coupling constant g in the Hamiltonian vanishes. We note that with 
£i{—Fz) = £^(Fz) and the symmetries of the auxiliary functions J-±{x,y) — ±T±{y,x), Siz becomes an odd function 
whereas Cxx is an even function of Fz- If eq. (|31|l was to hold, Dxx would have to vanish and Cxx would have to be 
a constant. This is true if and only if g — 0. Figs. [5] and O depict the dependences of Sz2 and of the correlation 
functions, respectively, on Szi for finite coupling strengths. In all cases, except for g = 0, and for the correlation 
Czz J the deviations from linearity are strikingly obvious. We note, however, that for small values of the Bloch vector 
component Siz the component S2z and correlation Cxx are almost constant and the other correlations Czz and Cyy 
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FIG. 3: The non- vanishing correlation functions Cxx, Cyy, Czz as a function of the Bloch vector component Sz for fie — 1 and 
(3g = 1.5. The correlation functions Cxx and Cyy clearly deviate from straight lines following from eq. For the correlation 

functions also an approximately linear regime exists in a neighborhood of Siz = 0. 

are linear. This behavior is in accordance with a convex blow-up map. Below we will come back to the blow-up in 
the linear response regime when the external fields are small. 

IV. DISCUSSION, IMPLICATIONS AND CONCLUSIONS 

We have illustrated Pechukas' verdict on the linear time evolution of open quantum systems [Tol | by a simple 
example. Moreover, we have demonstrated that the theorem holds for affine time evolutions: If two complete sets 
of pure system states can be prepared this more general class of evolution implies a factorizing preparation of the 
total initial density matrix where the environment density matrix must be independent of the system density matrix. 
Actually, Pechukas' original proof applies as well in the affine case. Nowhere in the proof he made explicit use of 
the homogeneity condition, i.e. R{Xps) ~ XR{ps), A real, that would render an affine R a linear map. 

At first, this may only seem a modest generalization of the original conclusion but it sheds some light on the crucial 
role of the inhomogeneity term of (generalized) master equations which appears when the initial density matrix does 
not factorize. The present analysis excludes that this term is independent of the initial reduced density matrix and 
actually is not merely an inhomogeneity of the otherwise linear master equation but must depend on the reduced 
density matrix in a nonlinear way. In those cases, when a Markovian dynamics is approached for long times this 
nonlinear term must vanish for sufficiently large times. It will do so, however, in a characteristic manner that depends 
on the particular initial reduced density matrix. 

In the present example only one set of pure system states, the eigenstates of CTz, are preparable and yet the affinity of 
the blow-up map implies a factorizing preparation. We note that in general the condition on the number of preparable 
pure states of the Pechukas theorem cannot be relaxed. A counter example is provided by the work of Karrlein and 
Grabert who considered a harmonic oscillator coupled to a bath of harmonic oscillators with a non-factorizing 
thermal preparation. This preparation allows for pure position states and still leads to a linear master equation. 

Another example of a non-factorizing preparation that leads to a linear master equation results from the following 
preparation procedure: (i) start with a factorizing density matrix pso Pbo at a time t = — <o, to > 0, (ii) turn on the 
interaction between the system and the bath, and (iii) use the density matrix that has evolved at t = as the result of 
the preparation. We term this preparation the factorize- and-wait preparation. The blow-up map R{ps) then assumes 
the form: 

R{ps) = e-^"'°Gi^\ps)pB^e^-"'\ (32) 
where Gt is the linear propagator of the reduced density matrix for the factorizing preparation, i.e., 

ps(i) =Gt(ps(0)) =trBe-i^*Vs(0)psoe^^*°. (33) 

Here, the inverse of the propagator Gt is needed in order to infer the proper system part of the factorizing density at 
t = —to from the density matrix ps that is to be prepared at t = 0. In view of possible fast relaxation processes and 
the build-up of system-bath correlations [23,|2J,|23 after the interaction has been switched on, this inverse propagator 
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will not be defined on the total set of possible density matrices [3^ . Still, it is a linear operator and thus the blow-up 
map (|32l) also is linear. It is obvious, however, that no pure states of the system can be prepared at time t = in 
this way, i.e., that the factorize-and-wait preparation does not provide pure states. Therefore, the conditions for the 
Pechukas theorem are not met and the theorem is thus not in conflict with the linearity of the blow-up map of the 
entangled factorize-and-wait preparation. 

Finally, we discuss the equilibrium preparation in the limit of weak external forces. In the region close to thermal 
equilibrium Mori's generalized quantum Langevin equations provide a proper description of the time evolution 
of the set of system operators that couple to the external fields. In this case, the preparable density matrices of the 
total system follow from eq. in linear approximation in the external fields Fj and hence assume the linear response 
form l34l: 



i Jo 

where p° denotes the equilibrium density matrix of the total system in the absence of external fields and {X)o the 
average of the operator X with respect to the density matrix p". The corresponding density matrix of the reduced 
system pg = Trs then becomes 



Ps-Ps + PT. f'dxTMpy-^ 



(X,-(X,)o)(p°rF,. (35) 



We recall that the external fields act on the system, i.e. that the conjugate operators Xi are system operators. The 
expectation values of these operators with respect to the density matrix pg are linear functions of the external fields 
by construction and can be written as: 

{X, - {Xj)o) EE TVS {X, - {X,)o) pf = E XJ^F^, (36) 

i 

where the response matrix Xij is obtained by inserting eq. (|34|) into the middle term of eq. (|36|1 . 

X^, f dx Tr (X, - (X,)) (p")!-- {X, - {X,)o) {p^f. (37) 
Jo 

Because Xij is an invertible matrix we may solve eq. (|35|l for the external fields Fi 

j 

In this way the external fields are expressed in terms of the linear functional {Xj — (X, )o) of the system density matrix 
p^. Using this relation in eq. (|34|) we find for the blow-up map the affine form R{pg) = Lp^ + x where 

LP^ = /^E fdx{py-^X,-(X,)o)ip°rx-\,TTsiX,~{X,)o)p^s, 

X = (39) 

With eq. (jSJ one obtains for the time evolution of the reduced density matrix 

p|(i) = Ttb U{t)L{pl)U\t) + Trs p° (40) 

where we used that the total density matric is invariant under the full time evolution. Actually this is a linear 
equation in the reduced density matrix ps- In order to see this one puts t = in cq. H4U|) . to express the trace of 
p" over the environment in terms of the reduced density matrix ps'. Tr^ p^ = ps — Tr^ L{ps). Hence, for the Mori 
preparation the time evolution of the reduced density matrix is linear. This is also in agreement with the findings for 
the above discussed model, see Figs. |21 and 13 
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APPENDIX A: CONVEX MAPS ARE AFFINE 



We prove that a difFerentiable convex map M from a Banach space Bi into a Banach space B2 is also affine. The 
derivative of M{x) at a; G i?i is defined as the hnear map DM [x) : Bi ^ B2 that is tangential to M at x. For further 
mathematical details see e.g. Ref. j35| : 

\\M{x + h)-M{x)-DM{x)h\\2 _p 
ll^lli^o \\h\\i 

where || ■ i = 1,2 denotes the norm in the respective Banach space. The convexity of M then requires that 
its domain of definition 'D{M) is convex, i.e. that with each pair x, y of elements of 'D{M) also all convex linear 
combinations, Xx + (1 — \)y, with < A < 1, belong to 'D{M). Moreover, this property is conserved under the convex 
map M: 

M{\x + (1 - \)y) = XM{x) + (1 - X)M{y). (A2) 
Taking the derivative with respect to A one finds from l|A2|l : 

DM{\x + (1 - X)y) {x-y) = M{x) - M{y). (A3) 

For A = one obtains 

M{x) = DM{y)x + M{y) - DM{y)y. (A4) 

The first term on the right hand side is linear with respect to x and the second and third term are independent of x 
and hence M is affine. 



APPENDIX B: PECHUKA'S THEOREM 



Pechukas proved in Ref. [W^ that a preparation is factorizing if the following conditions are satisfied: (i) the 
corresponding blow-up map of the preparation is convex; (ii) two different complete sets of pure system states can be 
prepared in the considered preparation class. 

The proof consists of two steps. The first step is a consequence of the positivity of the blow-up map R{ps) and 
makes use of the possibility to prepare pure states of the system. Assume {\ip)} is a pure state of the system that can 
be prepared. In the first step of the proof it is shown that the density matrix p of the full system has the form 

p = i?(|V)(^|) = |^)(V'lx (Bl) 

where the bath-density matrix x in general depends on the system state {ip}- 

For the proof we note that any density matrix of the total sytem can be represented as a weighted sum of products 
of pure system density matrices \ipi){ipi \ and bath density matrices Xi, 

P = Y.C^\2P^){i:^\x^, (B2) 
i 

where the states are normalized and orthogonal on each other: (ipil^j) — Sij, and the coefficients c; are not 
negative and add up to unity: J^i'^-i = 1- Taking the trace of R{\^^:) (ipl) over the bath one recovers the pure state 
E^nd on the other hand, using Tr^ Xi = 1 one finds from cq. ||B2|) : 

i 

Hence, all but one coefficients Ci vanish and eq. (|B1|I holds. 

Up to this point we have not made use of the convexity of the blow-up map and therefore the reduced bath density 
matrix x may still depend on the system state {ip). 

In the second step of the proof we first take two pairs of orthonormal system states, \tpi), IV'2) and \(pi), \ip2) that 
span the same two-dimensional subspace: 

|V^l)(V^l| + \^2){^2\ = + \V2){'P2\ = P, (B4) 
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where P is the projection operator on this subspace. The absolute values of the mutual scalar products then are 
determined by an angle a: 

{ipi,ipi) f = cos^a, I {^1,^2) f = sin^a, 

(V'2,</'i) P = sin^a, I (V'2, V2) P = cos^a. (B5) 

We assume that the two pairs of states ipi and ipi can be prepared. According to the first step of the proof the full 
density matrix that corresponds to either of the pure states is a product of the pure state and a density matrix of the 
bath: 

R{m{^P^\) = \^^){4'^\X{4'^), R{\^r){^r\) = \ip^) {ip^\xiV^) , « = 1, 2, (B6) 

where x(V'i) a-nd xiVi) ^^re bath density matrices. We shall show that they all are the same. For this purpose we 
make use of the convexity of the blow up map and consider its action on the density matrix that is proportional to 
the projection P onto the two dimensional subspace, ps — \P ■ The action of the blow up map can be represented in 
either of two ways by convex combinations of pure states, see eq. (|lj4|) : 

IV'l)('0l|x(V'l) + |'02)(V'2|x('/'2) = |<^l)(<^l|x(<yfl) + |V52)(¥'2|X(<^2)- (B7) 

Calculating the matrix elements with the four pure states V'ij '^ii one obtains the following equations relating the 
reference density matrices xi'4'i) ^^^d xi'fi)'- 

( Xi^i) \ _ ( |cosa|2 Isinap \ / x(^i) \ /p,o^ 

X(V'2) )'\\ sinap I cosap ) \ x{V2) ) ' ^^^^ 

and 

( Xi^i) \ = ( |cosa|2 Isinap \ f x(V'i) \ 

\X{^2) ) l^lsinap \cosa\'^ ) \x{'4^2) ) ' 

We may exclude the trivial cases when the pairs and (pi coincide and either cos^ a = 1 or sin^ a — 1. In all other 
cases, the four equations only have a solution with x(V'i) = x(V'2) = x('/'i) — x('/'2)- Hence, the bath density matrix 
is identical for all preparable system density matrices in the considered subspace such that a factorizable state of the 
total system results on this subspace. 

In order to apply this argument to higher dimensional system Hilbert spaces it must be possible to prepare sufficiently 
many pure system states. Starting as above with a two dimensional subspace that can be spanned by two pairs of 
preparable bases {V'ijj {Vi}^ * = 1, 2 one next considers the subspace spanned by either j-f/'i) and jV^a) or by \Lpi) and 
1^93) and finds x(V'i) = x(V'3) for the reference bath operators in the density matrix of the full system. That means 
that two different sets of pure states spanning the system's Hilbert space must be preparable. 
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